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We wish to call attention to a result which, to our 
knowledge, has not previously been noticed: in a BCS 
(Bardeen, Cooper and Schrieffer) superconductor |l| a 
range of pair end states may be accessed via both nor- 
mal and umklapp scattering. The situation resembles a 
double-slit experiment in which electrons have two paths 
to access the screen and cancel each other periodically in 
the production of the analogue of Young's fringes. 

Umklapp scattering is compatible with the canonical 
transformation 2] which is the basis of the BCS theory 
Discussion about umklapp contributions to superconduc- 
tivity started almost immediately after the BCS theory 
was forged [E 3 . This discussion has gathered renewed 
momentum after the high Tr superconductors have been 
discovered [E IE IE S fioj . It is therefore timely to 
discuss our new result in some detail. 

We show that the trial electron wave function in the 
standard BCS theory is not designed to accommodate 
states with multiple access paths. We find that, if we al- 
low all the end states of both normal and umklapp scat- 
tering to be generated by this wave function, then we are 
not able to normalize it. 

We extend the approach of BCS, who already modified 
their T = trial function to accommodate single particles 
and excited pairs that appear at T > 0. We modify 
the BCS trial function at both T = and T > to 
include occupation of states by both normal and umklapp 
scattering. We show that the modified trial function is 
always normalized. 

We also show that the effect of umklapp scattering is to 
cancel the pairing effect of normal electron-phonon scat- 
tering. The cancellation is actually the physical founda- 
tion for the average phonon frequency significantly lower 
than the Debye frequency, introduced by BCS. We re- 
strict our discussion exclusively to the BCS formalism, 
based on the principle of variation. We do not discuss 
for example the Eliashberg formalism which derives 
from a different principle. 



In the BCS theory the standard procedure of variation 
is followed. It starts with the reduced Hamiltonian 



(1) 



k,k' 



where Vk,k' measures the strength of the electron-phonon 
interaction in second order, 



a -kl a kT 



(2) 



generates pairs, a + being the particle generation oper- 
ator, k electron wave vector, f and j spin. In Eq. JU 
k runs over each and every state of the first electron 
Brillouin zone exactly once with spin f (we assume a sin- 
gle band throughout our discussion). Simultaneously — k 
runs over the same range with spin {. An alternative ap- 
proach, which we do not adopt, is to let k in Eq. Q run 
over half the Brillouin zone and sum the spins. Either 
way, we may have spin j and J, for any specific value of k. 
Consequently the paired particles may be in (k | , — k J.) 
or(k|,-kt). 

In Eq. Q the electron term is in and b^, which is not 
entirely equivalent to the standard electron Hamiltonian 
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(3) 



The electron term in Eq. will vanish when applied to 
single particles, but the Hamiltonian in Expression © 
will not: Eq. Q is designed only for a pair ensemble. In- 
deed even Expression J2J is designed only for electrons in 
crystals, where we have discrete electron energy ek, which 
is produced automatically by the Hamiltonian when an 
electron is detected by the particle destruction operator 



III. TRIAL FUNCTION 

Next is introduced the trial wave function for quasi- 
particles. At T = this function is of the following form 
in the BCS theory Q: 



(4) 
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/ik is the pair occupation probability and |0) the pair 
vacuum. Here we seek the equilibrium configuration of 
the particle ensemble under the influence of the electron- 
phonon interaction. The trial function must accommo- 
date all the initial and end states of the paired quasi- 
particles, which are indicated by k and k' in Eq. with 
occupancy probabilities to be determined for both the 
initial and end states in a single step of the process of 
variation. 

Now in Eq. the product over k also runs over the 
first electron Brillouin zone, i.e. the number of orbits in 
Eq. (0J equals the number of Bloch states. However, 
this does not mean that we cannot add further orbits to 
Eq. (@J . The growth of the superconductive energy gap 
can be compared with the splitting of the spectral line 
of say the hydrogen atom, where the electron has several 
configurations, with different patterns of particle loca- 
tion probability and/or spin but the same energy. This 
energy degeneracy splits to reveal further structure in the 
presence of e.g. an external magnetic field. Similarly, a 
particle pair has two configurations, generated by 

vi-^k+y^^ (5) 

and 

Vl-^k b+ - y/k^ (6) 

respectively, whose energy splits from the Bloch energy ek 
in the presence of the electron-phonon interaction. The 
energy of the pair generated by Expression (J5J shifts 
downwards to give the energy gap. The energy of the 
pair generated by Expression © shifts upwards. Appar- 
ently, in addition to the Bloch orbits, we may have two 
more orbits per primitive cell. Eq. does not include 
all the available orbits, because the paired states there 
suffice to accommodate the physics: the superconductive 
current is carried only by the paired quasi-particles with 
energy shifted downwards under the electron-phonon in- 
teraction. 

When T > 0, Eq. (@J no longer suffices to accom- 
modate the physics. Indeed, if we continue to employ 
Eq. J3J as our trial function, then we have no chance to 
see what happens when superconductive pairs break into 
single particles or are excited under the attack of ther- 
mal phonons. To cope, BCS modified Eq. Q into the 
following form: 

Wexc = n 4-. n (v/^kT b+ 

k"(5) k'CP) 

-v^) n (vi-/ik+ v^k"&k)io> (7) 

HS) 

where S, V and Q specify the states occupied by single 
particles, excited pairs and ground pairs respectively, a 
(could be denoted as Ok) stands for either \ or \ depend- 
ing the value of k 0]. It is interesting to note that in 
Eq. 10 k, k' and k" do not have to be different. In- 
deed we can check through direct calculation that, when 




FIG. 1: A spherical Fermi surface (open circle) and a pair of 
electrons in states k and — k. This pair is scattered into states 
k + q + G and — k — q — G in the umklapp process, q being 
the phonon wave vector, G the reciprocal lattice vector. The 
shaded circles represent two identical spherical phonon zones. 
We use angles a and f3 to measure the size of the sections 
of the Fermi surface being intersected by these phonon zones. 
We have a > /3 because a larger section of the Fermi surface is 
intersected by the phonon sphere centered at this surface (the 
other sphere is centered above the Fermi surface at height go). 
As a result, the electron state — k q G must lie inside the 
phonon sphere centered at k, where all the states are involved 
in normal scattering. When k + q + G runs over the shaded 
phonon zone on the left, — k q G also runs over a sphere 
indicated in part by the broken white circle (centered above 
the Fermi surface at height go to mirror the phonon sphere 
on the left). 



k = k' = k", the wave functions of the single particle, 
the excited pair and the ground pair are orthogonal to 
each other: the total number of states in Eq. (J7J) can be 
three times as large as that in Eq. J2J. 



IV. DESTINATION STATE WITH MULTIPLE 
ACCESS PATHS 

Consider in Fig. ^ a spherical Fermi surface that does 
not intersect the boundary of a Brillouin zone: suppose 
also there is a minimum separation qo between adjacent 
Fermi surfaces. An electron in state k may be scattered 
into k + q to lie on the adjacent Fermi surface if the 
phonon momentum q exceeds go ■ Since it makes no phys- 
ical difference if the reciprocal lattice vector G (parallel 
to k for simplicity) is added to an electron state, we can 
replace k + q with k + q + G, which lies on the original 
Fermi surface. This state is known as the end state of 
umklapp scattering and, in a non-superconducting metal, 
cannot be accessed through normal scattering (in the re- 
duced zone scheme), at least for a range of q values. 

In a superconductor one destination state may have 
multiple access paths. For example in Fig. ^k and k 
are paired and, correspondingly, k + q + G and k — 
q — G are also paired: scattering of the particle pair is 
synchronized, reminiscent of the synchronized motion of 
rotating doors. In the umklapp process in Fig. ^ we have 
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the following possible chains of events concerning state 
occupation: 



k+q^k+q+G^kqG 



(8) 



The first two arrows represent scattering and umklapp 
folding. The third arrow denotes synchronized occupa- 
tion of the state — k — q — G by the particle initially in 
state — k. These events can be compared with say pas- 
sage of an electron through a slit, or reflection of a photon 
by a mirror, in an interference experiment. It is evident 
from Fig.Q]that — k — q— G lies inside the phonon sphere 
centered at k, where all the electron normal-scattering 
end states lie, and therefore we also have the following 
normal-scattering event: 



G 



(9) 



Apparently the state — k q G can be accessed via two 
different paths. 

We use Af and U to denote collections of all the (sin- 
gle) end states of normal and umklapp scattering, respec- 
tively, which form thin layers about the Fermi surface, in 
which empty states are available to accommodate scat- 
tered particles. In Fig. [2] Af is a section of the Fermi 
surface intersected by the phonon sphere centered at k, 
whose radius is qo (Debye momentum) to include all the 
phonons allowed in the Debye model. This phonon sphere 
may be transported to the left by a lattice constant G. 
The transported phonon sphere also intersects a section, 
U, of the Fermi surface to accommodate particles scat- 
tered from k in the umklapp process. We separate Af 
into two sub-sets, Afi and A2, where Afi accommodates 
particles scattered from — k in the umklapp process. It 
is apparent from Fig. [21 that an extensive range of states 
in Af can be accessed via both normal and umklapp scat- 
tering. 

We also use Af and U to specify paired end states of 
normal and umklapp scattering (or normal and umklapp 
pairs for short). In Fig. [5J the single particle in state 
k may have spin f or J,. Consequently the normal or 
umklapp pair may be in (k' f, — k' J.) or (k' |,k' f) after 
scattering. The shape of Af and U is the same in both 
cases. 



V. LIMITATION OF BCS TRIAL FUNCTION 

The BCS trial function is not designed to accommo- 
date states with multiple access paths, such as those de- 
scribed in Expressions (JSJ and ©. To appreciate this, 
we study the example in Fig. [21 where we illustrate how 
states are linked by Vk,k' in Eq. when k is in the di- 
rection of G. We write Eq. (0J in the following equivalent 
form 




FIG. 2: A spherical Fermi surface (open circle) with two 
phonon spheres of radius qo (Debye momentum). In umklapp 
scattering the particle in k is scattered into U, which is a 
layer about the Fermi surface shown as a thick black arc. In 
normal scattering the particle in state k is scattered into Af = 
Afi + A/2 , which is also a layer about the Fermi surface shown 
as a thick arc in black and white, Afi and U are symmetric 
with respect to the center of the Fermi sphere. The sizes of Af 
and U are measured by A and B, which are maximum values 
of a and (3 in Fig. Q There is a hole in the phonon spheres 
to indicate that umklapp scattering is not invoked when the 
phonon momentum is less than qo. 



with 



1*)= n n n (v^~~^+ 

k'(,A/i) k'(jV 2 ) k'(W) 



h k , 6+)|$) (10) 



k"(7J) 

where 7Z specifies states not in Afi, Afi and U. By defini- 
tion for each k' in IA we have — k' in Afi and vice versa, 
so that we can write the product over k'(A/i) in Eq. ltTU|) 

as 

n (vi-zi-k'+v^frtk') (u) 

Combining Eqs. (|10|l and (|11|) . we find 

(*i#)= n [1+2VM1-M 

V(U) 

xV/i-k'U-ft-k') (12) 

so that > 1 even when ($|$) = 1. Obviously 

the standard BCS theory is not designed to endorse the 
umklapp process in Fig. ^ 

We cannot attribute the self-contradictory result in 
Eq. (fP2|> to counting k + q and k + q + G in Fig. ^ as 
different states, because the states in Eq. Ill 01) are all in 
the first electron Brillouin zone (Fig. El. We might like 
to attribute Eq. lfP2|) to the proposition that normal and 
umklapp scattering drive two pairs into the same state, 
and thus violate the exclusion principle. However, the en- 
tity scattered in both the normal and umklapp processes 
is the same pair in k and — k (Figs. ^ El : we do not 
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have two pairs to violate that principle. It appears that, 
in Eq. (0J, the number of orbits is simply not enough 
to cope with the situation that some areas of the Fermi 
surface, such as Af\ in Fig. |3 suffer enhanced high prob- 
ability of pair bombardment when normal and umklapp 
scattering coexist. Indeed Eq. (|12(l does suggest that the 
number of pairs landing in U and Af\ exceeds the limit of 
its design. 



VI. MODIFIED TRIAL FUNCTION 

We extend the approach of BCS, who modified the 
trial function at T = 0, Eq. into the trial function 
at T > 0, Eq. @, to accommodate single particles and 
excited pairs. We will modify the trial function at both 
T = and T > to accommodate umklapp scattering. 
We will add more orbits to Eq. to allow umklapp 
pairs to land in say U and Af\ in Fig. |2 without self- 
contradiction as in Eq. I|12|). Since Expression JSJ has 
already been used to generate normal pairs in Eq. 
we have no choice but to use Expression © to generate 
umklapp pairs and find 



k(U) 



m n (vi-^+v^^io) (13) 

k(AQ 

as our trial electron wave function in the presence of umk- 
lapp scattering. 

According to Eq. (fC^I the equilibrium state of the par- 
ticle ensemble consists of a mixture of normal and umk- 
lapp pairs. We find through direct calculation that we 
always have | = 1 if we can prove 

hk = h-u (14) 

where k and — k are specified by IA and A/", respectively, 
in order to avoid a self-contradiction similar to that in 
Eq. 1|12[) . Eq. (|14fl must be true, because both the Fermi 
sea and reciprocal lattice are symmetric with respect to 
k = 0: interchanging k and — k makes no difference in 
Eq. the Hamiltonian. 

We have a more serious problem: in Eq. i|13|l we have 
actually assumed that /i k depends only on the numerical 
value of k, regardless of whether k is specified by A/" or 
hi. Since the condition to find ft, k may be different in 
Af and U, it is still possible that Eq. (|13|) could be self- 
contradictory, as was Eq. Q , in the presence of umklapp 
scattering. We will address this in the next section. 



VII. OCCUPANCY OF NORMAL AND 
UMKLAPP PAIRS 

We use Eq. (|13fl to evaluate the energy of the pair 
ensemble 



where Hbcs is defined in Eq. QJ. We find 
W = 2 £ k^k + 2 J2 £ k(l - fck) 

k(AT) k(U) 
-W n ~ W 22 + Wi2 + W 2 1 



(16) 



with 



Wi 2 



J2 E ^WMi-^OVMi-M 

k(Af) k'(A0 

Vk.kVMi - fckVM 1 - foe) 

k(jV) k'(W) 



W21 and W22 are identical to W12 and Wn, respectively, 
save that Af is replaced by U and vice versa. 

We minimize W in Eq. fTS} with respect to h^. When 
k is specified by Af we have 



dW 



2e k 



1-2/lL 



^/i k (l - /ik) 



E-E 

k'(AT) k'(U) 



Vk,kV^k-(l-/ikO (17) 



Here and are linked together by the constraint 



hu = - 



£k 



E(k) 



(18) 



with E(k) = v / A 2 (k) + el. Letting dW/dh k = 0, we 
find from Eqs. (|T7|> and (|T%|) 



A(k) 



E-E 

k'(AT) k'(U) 



Vic,] 



A(kQ 
2£(k') 



(19) 



as our self-consistent equation, which is identical to the 
standard BCS self-consistent equation 0, save the terms 
in k'(U), which reduce and may even cancel the effect of 
the terms in k'(Af). 

When k is specified by U, we have 



dW 



-2e k 



1 - 2fr k 
y/hk(l - /i k ) 



E-E 

k'(JV) k'(U) 



V Kk .y/h k ,(l-h w ) (20) 



W = (^IHbcsI*) 



(15) 



which also leads to Eq. (|19|l . Apparently we have the 
same A(k) for the same numerical value of k, regardless 
of whether k is specified by Af or U. This A(k) leads 
through Eq. I|18|) to ft. k , which also depends only on the 
numerical value of k, and thus ensures 1 = 1 for 
Eq. 1131). 



VIII. FRUSTRATION OF 
SUPERCONDUCTIVITY 

The principle of variation can be compared with shak- 
ing, which shapes a physical system into a configuration 
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(a) 



(b) 



(c) 



FIG. 3: Configurations of (a) bar with minimum energy, (b) 
bar with maximum energy and (c) two bar system. 



with minimum energy, such as the horizontal orientation 
of the bar in Fig. [3] (a) in a gravitational field. How- 
ever the process of variation also produces a configuration 
with maximum energy, such as the vertical orientation of 
the bar in Fig. [31 (b) . Although we should choose the 
configuration with minimum energy whenever possible, 
some physical systems may frustrate us from doing so. 
For example, when the two bars are joined together to 
become an 'L' in Fig. [21(c), we cannot minimize their en- 
ergy as readily as we do in Fig. [21(a), because minimum 
energy of one bar implies maximum energy of the other, 
and vice versa. 

We have a similar situation when normal and umk- 
lapp scattering coexist. The pair occupation probability, 
/ik, has a phase-like property, which becomes clear when 
we let /ik = cos 2 #k, so that the normal pair generation 
operator in Expression JSJ) becomes 



sin 6>k + cos 



(21) 



The resultant wave function is normalized at any value 
of #k- However, in order to maximize the energy gap, we 
have to let 



cos 2 9i, 



1 



1 



(22) 



which is equivalent to Eq. (|18|l . We can also write the 
umklapp pair generation operator in Expression © as 



sin(0 k - 90°) + cos((9 k - 90°) 



(23) 



It is easy to check that all the results here or in the 
BCS theory will not change if we replace the —90° phase 
difference in Eq. |23J> with 90°. 

According to Eq. I|19|l the normal pair is in a phase 
to sustain superconductivity, but the umklapp pair is in 
a phase to frustrate it. However, if we replace Eq. i|18|) 
with 



h v = - 



which can be written as 



1 + 



^k 



E(k) 



cos 2 (0 k -9O°) = ~ 



1 - 



£k 



(24) 



(25) 




FIG. 4: A spherical fermi surface (open circle) and two identi- 
cal phonon spheres, G being the lattice constant. Phonons in 
the unshaded circles do not invoke umklapp scattering. The 
particle in state k is scattered into k + p and k + q + G in 
the normal and umklapp processes, respectively, by phonons 
with the same momentum. Energy of the end state in both 
processes is at the same level above the Fermi surface (bro- 
ken arc, with an exaggerated constant height above the open 
circle) . 



then in Eq. I|19|) Af will be replaced by U and vice 
versa: umklapp pairs sustain superconductivity but nor- 
mal pairs frustrate it. Apparently both normal and umk- 
lapp pairs alone may lead to an energy gap, but their 
coexistence mutually cancels their effect. 

The 90° phase difference in Eq. I|22l) . reminiscent of 
the angle between the two bars in Fig. 0H(c), reveals an 
interesting dilemma. On the one hand the normal and 
umklapp pairs have to be locked 90° apart form each 
other in phase, otherwise our trial function, Eq. Q13p. 
cannot be normalized. On the other hand, we have to pay 
the price that either normal or umklapp scattering will 
be in a phase to frustrate superconductivity. Apparently 
only a fraction of phonons may survive cancellation to 
play their traditional role in the standard BCS theory 
and this, we suspect, was actually implied by BCS some 
50 years ago. 



IX. AVERAGE PHONON FREQUENCY 

In Eq. (|19fl there can be no term in k'(W) to survive 
cancellation by a term in k'(A/"), at least for a spherical 
Fermi surface that does not intersect the boundary of a 
Brillouin zone, and with the constraint in Eq. (|18|) . To 
see this, let k be a state on the Fermi surface (not in the 
direction of G for generality, see Fig. 0} . We use k + p 
to denote the end state of normal scattering, which is 
confined within the shaded phonon sphere, of radius p 
and centered at k, as a result of momentum conserva- 
tion. On the other hand, energy consideration requires 
that k + p must be at a level near the Fermi surface, 
which is marked as a broken arc in Fig. 01 somewhat ex- 
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aggerated in its height above the Fermi sphere. There is 
an unshaded hole inside the phonon sphere (radius ^ q Q , 
because k and G are not parallel) where umklapp scat- 
tering is not invoked. We assume that k + p is not in the 
hole. 

If we transport the phonon sphere to the left by a re- 
ciprocal lattice constant G, then we find all the end states 
of umklapp scattering allowed by momentum conserva- 
tion (the hole is completely outside the Fermi sphere but 
contacts its surface). ^From the intersection between the 
surface of the transported phonon sphere and the broken 
arc in Fig. 0] we find a state k + q + G with 



q = p, 



Ck+q+G — £k+p 



(26) 



i.e. both the phonon momentum and end state energy 
are the same in both normal and umklapp scattering. In 
Fig. 0] we have k = (k x , k y ,0). When k z ^ we can still 
find k + q + G to satisfy Eq. I|26|l by moving the phonon 
sphere in a proper direction over a proper distance. 

We are reminded that, starting from p, we can find 
more than one q to satisfy Eq. I|26|) . However, we do have 
a one-to-one relation between p and q. To prove this, we 
reverse the process in the previous two paragraphs. We 
start from any end state of umklapp scattering, k+q+G, 
and then find k + p to satisfy Eq. Apparently p 

will not fall into the unshaded hole in the phonon sphere: 
this one-to-one relation is between the states in U and 
those states in Af which are subject to both normal and 
umklapp scattering. 

According to Eq. (|26[1 we have 



Vk, k+q+G = Vk,k+p 



where p| 



Vk,k+p 



2hw p \M p \ 



(H) 2 ~ ( £ k+p - £ k) 2 



(27) 



(28) 



M.p being the matrix element, lo p phonon frequency, 
both functions of p (Debye model), and Vk,k+ q +G is de- 
fined similarly. For an isotropic energy gap we also have 



A(k + q + G) = A(k + p) 



(29) 



and hence E(k + q + G) = E(k + p). Eqs. J23 and O 
lead through Eq. (JT^ to 



A(k) = V W 



k'(B) 



A(k') 
2E(k') 



(30) 



B specifies states that do not invoke umklapp scattering, 
which is confined within the unshaded circle around k in 

Fig. a 

We should remember that the Fermi sea and the recip- 
rocal lattice have incompatible symmetries. Therefore in 
Eq. I|30|) B and hence A and E depend on the direction 
of k even when the Fermi sea is spherical. For simplicity, 



we have ignored the k-dependence of B, as a basis for 
deriving Eq. (|29|l from Eq. i|26|) . In future work, we will 
have to generalize our proof in the previous paragraphs 
by replacing Eqs. |J27J| and with 

A(k + q + G) A(k + p) 

V^k+q+G-— ^— — — - Kk,k+p OEV1 , | ^ K 61 ) 



'2P(k + q + G) 



2P(k + p) 



as our condition to seek q for a given state p (and vice 
versa). We will have to find an approximate value of q 
through the procedure in Fig. Then we will adjust q 
slightly (assuming that phonon states are sufficiently nu- 
merous for small adjustment) to balance Eq. (|31|l . which 
also leads to Eq. . 

Eq. (|30|l is in accord with the view that the coexistence 
of normal and umklapp scattering frustrates supercon- 
ductivity. Indeed, if we cut the phonon frequency off at 
the threshold of umklapp scattering, and go through the 
standard BCS formalism, then we find a self-consistent 
equation in exactly the same form as Eq. (|30|l . 

BCS introduced ui, the so-called 'average phonon fre- 
quency', to define the range of phonons when they in- 
tegrated the self-consistent equation pj. In the recent 
literature a liberty is often taken to replace u> with the 
Debye frequency lod- However we see from Eq. <|30ll that 
the BCS average phonon frequency is based on a very 
deep insight. Numerically w is likely to be significantly 
lower than up. According to Figs. and the phonon 
momentum in Eq. 13U|) cannot exceed q$ when k is in the 
direction of G. For a bcc lattice with valency 1 we have 
qo/^F = 0.267, kp being the radius of the Fermi sphere, 
or qo/qD = 0.212: when k is in the direction of G the 
average phonon frequency of BCS is merely ~ 20% of the 
Debye frequency. 



X. EXCITED STATE 

For the convenience of the reader, we explain the BCS 
theory at finite temperatures briefly and then extend the 
formulation to include umklapp scattering. BCS divided 
reciprocal space into cells, Ak, containing N k pair states. 
There are Sk single particles and Pk excited pairs in Ak, 
other states being occupied by ground pairs, giving 



N u -S k - P k 



(32) 



as fractions of Ak occupied by single particles, excited 
and ground pairs, respectively. The kinetic energy of the 
particle ensemble, found from Eq. J7J, the electron term 
of Eq. QJ and the Bloch energy contribution, is 



W K e = J! £k 



-2hi 



Nu~Su- Pk 
A^k 



(33) 
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^KE=2^e k [/ k + (l-2/ k )fc k 



with 



/k 



S k P k 

2N k 7V k 



(34) 



(35) 



In Eqs. (|33|l and Il34() e k is averaged in Ak (a cell), k 
runs over the cells (one k for each cell). The interaction 
energy is found to be 



Wi 



J2 wmi-^p - 2/ k ) 



k,k' 



xVftk'(l-hk')(l-2/k') 
The entropy is given by 



(36) 



5 = -2fc B ^ [/ k In / k + (1 - / k ) ln(l - / k )] (37) 

k 

ks being the Boltzmann constant, which leads to the free 
energy 



F = W KE + W l -TS 



(38) 



We have 
OF 



x X!Vv / M^ T M(l-2/k-) (39) 



which can be compared with Eq. I|17|l . When dF/dh^ = 
0, Eq. jSHJl leads through Eq. (JTHJ to 

A(k)=|:^^^(l-2/ k (40) 

It is interesting that k' runs over the same range at any T. 
Therefore, when the temperature rises, ground pairs are 
replaced by single particles and excited pairs gradually 
in each and every cell, Ak', rather than being driven out 
altogether cell by cell. We also have 



dF f 

— = 2E(k) + 2k B T\u 

djk 1 - Jk 



or 



1 — 2/ k = tanh 



2k B T 



(41) 



(42) 



when dF/dfk = 0. Combining Eqs. gUJ and g2J>, we 
find 

A(k).W kk ,^Uanh^ (43) 



which is the self-consistent equation at T > in the BCS 
theory [|. 

Eq. 1)43(1 . though sophisticated and elegant, manifests 
a simple idea, viz all states fall into two categories: those 
allowed at T = and those not, to evaluate the entropy 
in Eq. I|37|l . Once we endorse this idea, the rest of the 
derivation is automatic. For each k we have three con- 
figurations, whose contributions to the ensemble energy 
are well-defined. Thermodynamics governs the number 
of the particles in each of the configurations, which in 
turn governs the magnitude of the energy gap. Appar- 
ently, this simple idea will not change when umklapp 
scattering is included. 

In the presence of umklapp scattering, we replace 
Eq. Q with: 



l*>« 



n n <■>* n n 



k"(5i)k"(5 2 ) 



k'(Pi)k'(S 2 ) 



M II II Wl-h* + Vh^b+)\0) (44) 

k(Si)k(-p 2 ) 

where Q and V specify ground and excited pairs, respec- 
tively, S specifies single particles, as in Eq. J7J), indices 
1 and 2 specify the normal and umklapp processes, i.e. 
we separate N and U into Si, Vi, Qi and S2 , P2, G2 re- 
spectively. Note that, when ground pairs are generated 
by Expression ©, we use Expression @ to generate ex- 
cited pairs, similar to the arrangement in Eq. Q. On the 
other hand, when ground pairs are generated by Expres- 
sion ©, we use Expression © to generate excited pairs, 
which raises the energy of the pair ensemble, in recogni- 
tion that umklapp pairs make less negative contribution 
to that energy the fewer their number. 

We divide M and U into cells Aki and Ak2 respec- 
tively. We have 



fflk Piw ^Vlk — Sik — Pi k 
AW AW Nik 



(45) 



as fractions of Aki occupied by single particles, excited 
and ground pairs. We also have 

S2U P2W N 2 k — S 2 k — P 2 k 

AW AW A^2k 1 ' 

as fractions of Ak2 occupied by single particles, excited 
and ground pairs. Letting 



h = 



Pn 



S 2 k 



P2k 

N 2k 



(47) 



2JV lk 7V lk 2iV 2k 
for k specified by J\f or U respectively, we find 

W^KE - 2e k[/k + (1-2/,,)^] 

k(jV) 

+ 5Z2e k [/ k + (l-2/ k )(l-/j k )] (48) 
k(W) 

as the kinetic energy, which can be compared with 
Eq. (EH). We also find 



Wi = -Wn - W 22 + Wu + W 2 i 



(49) 
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with 



w u = E VA(i-MA'(i-V) 

k(AO k'(AQ 

x(l-2/ k )(l-2/ k ,) 



and 



w 12 = 5^ J] vi, k VMi-/ik)VMi-M 

k(JV") k'(W) 

x(l-2/ k )(l-2/ k 

W21 and W^22 are identical to W12 and Wu, respectively, 
save that Af is replaced by U and vice versa. 

Now we minimize the free energy F, which is defined 
by Eqs. fSJl, igEJ and (gSJ when k is specified by TV. We 
find 



dF 
dhir 



2e k (l - 2/ k ) 



1 -2/lL 



= (l-2/ k ) 



E-E 

k'(A0 k'(W)- 



Vk,kV^k'(l-/ikO (1-2/kO (50) 



which can be compared with Eq. i|39|l . When dF/dh^ 
0, Eq. (jSni leads through Eq. (JTSJ) to 



A(k) 



E-E 

L k'(A0 k'(uy 



which can be compared with Eq. Eqs. l)4*T)) and l|l2^ 

do not change in the presence of umklapp scattering, and 
this leads through Eq. I|51|) to 



A(k) 



E^E 

It' (AO k'(W). 



Jr A(k') , B(k') 
Vk k' - -\. 1 tanh ■ 



2£(k') 



2fc B T 



(52) 



as the self-consistent equation for T > in the presence 
of umklapp scattering. Eq. I|52f) can be compared with 
Eq. g5J. 

When k is specified by U we define the free energy as 



F 



-W KE -Wi-TS 



(53) 



which can be compared with Eq. 1 155) 1. Wke and Wi are 
defined in Eqs. I|48|l and l|49|l respectively, S is defined 
in Eq. fgft. In Eq. (£5) W = W K E + Wi is negative, 
because umklapp pairs make negative contributions to 
the energy of the pair ensemble. This means that we have 
SQ = —dW if the change in Q is caused by the change 
in the number and/or occupancy of umklapp pairs, Q 
stands for heat. Therefor we have to let F — —W — TS, 
in order to have dF = SQ-TdS-SdT = -SdT, which is 
the correct expression of dF (for constant volume) . Now 
Eq. lf5U)l becomes 



dF 



-2e k (l - 2/ k 



1 - 2hv 



= (l-2/ k ) 



E-E 

k'(A0 k'(W)' 



Vk,kVMl-M (1-2/kO (54) 



which also leads to Eq. lfST)l . We still have Eqs. l(TT)) and 
g2Jl, which lead through Eq. (£>TJ to Eq. . 

If we adopt the view that superconductivity is frus- 
trated when normal and umklapp scattering coexist, then 
we can go through the standard BCS formalism at T > 
to find 



A(k) 



k'(B) 



A(k') , E{k') 
k' - =r. - tanh - v 



2S(k') 



2fc B T 



(55) 



which is identical to Eq. (|52|l . S (defined in Section Hx)l 
specifies the range of k', which does not depend on T, 
similar to the range of k' in Eq. 14U|) or (|43|l . 



XI. CONCLUSIONS 

In the conventional picture of single electron scattering 
in a metal the normal and umklapp scattering processes 
are associated with distinctively different ranges of mo- 
mentum transfer. There is no possibility for these two 
processes to interfere with each other in the sense that, 
with the same initial state, the end states of normal and 
umklapp scattering do not coincide. In this picture umk- 
lapp scattering is effectively treated as normal scattering 
in an extended range of momentum transfer and, when 
applied in the BCS theory, it tells us that umklapp scat- 
tering also lowers the ensemble energy and thus enhances 
superconductivity 

We recognize that, when the scattered entity is not a 
single electron but a pair, the end states of normal and 
umklapp scattering may coincide and the two processes 
may interfere. We find a formal deficiency in the stan- 
dard BCS ground state wave function, which is designed 
to count the statistics of scattering of pairs in the normal 
or umklapp process alone and leads to a self-contradiction 
when these two processes coexist. A suitable trial func- 
tion locks the end states of normal and umklapp scat- 
tering 90° apart in phase but leads to the consequence 
that superconductivity is frustrated when the normal and 
umklapp processes coexist. Then only a fraction of the 
phonons actually takes part in electron pairing. Our con- 
clusion is reached within the BCS formalism which is 
based on the variational principle. 

Another approach, the Eliashberg formalism, is based 
on the one-electron Green's function, and here too su- 
perconductivity is enhanced when umklapp scattering is 
introduced [13) . However, in this formalism, which is 
non- variational and does not involve a trial function, the 
physics of pair scattering and interference effects may be 
manifest in a manner different from that in the BCS vari- 
ational approach. Further consideration of the contribu- 
tion of umklapp phonons in different formalisms would 
appear to be merited. 
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